Tests of relativity using a microwave resonator 
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Abstract 

' The frequencies of a cryogenic sapphire oscillator and a hydrogen 

, maser are compared to set new constraints on a possible violation of 

Lorentz invariance. We determine the variation of the oscillator frequency 
as a function of its orientation (Michelson-Morley test) and of its velocity 
J£T , (Kennedy- Thorndike test) with respect to a preferred frame candidate. 

I ■ We constrain the corresponding parameters of the Mansouri and Sexl test 

theory to 5-/3+1/2 = (1.5±4.2) x 1(T 9 and/3-a-l = (-3.1±6.9) x 1(T 7 
f"*) , which is equivalent to the best previous result for the former and repre- 

■ sents a 30 fold improvement for the latter. 

PACS number(s): 03.30.+P, 06.30.Ft 

q-i The Einstein equivalence principle (EEP) is at the heart of special and gen- 

I ■ eral relativity [1] and a cornerstone of modern physics. One of the constituent 

£jf)! elements of EEP is Local Lorentz invariance (LLI) which, loosely stated, pos- 

• • ■ tulates that the outcome of any local test experiment is independent of the 

. £^ | velocity of the (freely falling) apparatus. The central importance of this postu- 

■ late in modern physics has motivated tremendous work to experimentally test 
LLI [1]. Additionally, nearly all unification theories (in particular string theory) 
violate the EEP at some level [2] which further motivates experimental searches 
for such violations of the universality of free fall [3] and of Lorentz invariance 
[4, 5]. 

The vast majority of modern experiments that test LLI rely essentially on 
the stability of atomic clocks and macroscopic resonators [6, 7, 8, 9], therefore 
improvements in oscillator technology have gone hand in hand with improved 
tests of LLI. Our experiment is no exception, the 30 fold improvement being a 
direct result of the excellent stability of our cryogenic sapphire oscillator. Addi- 
tionally its operation at a microwave frequency allows a direct comparison to a 
hydrogen maser which provides a highly stable and reliable reference frequency. 

Numerous test theories that allow the modeling and interpretation of exper- 
iments that test LLI have been developed. Kinematical frameworks [10, 11] pos- 
tulate a simple parametrisation of the Lorentz transformations with experiments 
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setting limits on the deviation of those parameters from their special relativistic 
values. A more fundamental approach is offered by theories that parametrise the 
coupling between gravitational and non-gravitational fields (THe/i [12, 1, 13] or 
Xg [14] formalisms) which allow the comparison of experiments that test differ- 
ent aspects of the EEP. Finally, formalisms based on string theory [2, 3, 4] have 
the advantage of being well motivated by theories of physics that are at present 
the only candidates for a unification of gravity and the other fundamental forces 
of nature. 

Owing to their simplicity the kincmatical frameworks of [10, 11] have been 
widely used to model and interpret many previous experiments testing LLI [6, 8, 
9, 19, 20]. In order to compare our results to those experiments we will follow this 
route in the present work (an analysis of our experiment in the light of other test 
theories being relegated to a future publication). Those frameworks postulate 
generalized transformations between a preferred frame candidate E(T, X) and 
a moving frame S(t, x) where it is assumed that in both frames coordinates are 
realized by identical standards (e.g. hydrogen masers for the time coordinates). 
We start from the transformations of [11] (in differential form) for the case where 
the velocity of S as measured in S is along the positive X-axis, and assuming 
Einstein synchronization in S (we will be concerned with signal travel times 
around closed loops so the choice of synchronization convention can play no 
role): 

™ 1 / , vdx\ dx v ( , vdx\ „„ dy ,„ dz 

d T=-(dt + ^);dX = - + -(dt + ^);dY=f,dZ=- (1) 

with c the velocity of light in vacuum in S. Using the usual expansion of 
the three parametrs (a w 1 + av 2 /c 2 + £>(c~ 4 );6 « 1 + f3v 2 /c 2 + £>(c~ 4 );d « 
1 + Sv 2 /c 2 + C(c- 4 )), setting c 2 dT 2 = dX 2 + dY 2 + dZ 2 in S, and transforming 
according to (1) we find the coordinate travel time of a light signal in S: 

dt = ^ ( 1 - (/? - a - 1) ^ - Q - + S\ sm 2 9^pj + 0( C - 4 ) (2) 



where dl = \J dx 2 + dy 2 + dz 2 and 9 is the angle between the direction of light 
propagation and the velocity v of S in S. In special relativity a — —1/2; (3 = 
1/2; 6 — and (1) reduces to the usual Lorcntz transformations. Generally, 
the best candidate for S is taken to be the frame of the cosmic microwave 
background (CMB) [15, 16] with the velocity of the solar system in that frame 
taken as v Q w 377 km/s, deck w -6.4°, RA w 11. 2h. 

Michclson-Morley type experiments [17, 6] determine the coefficient Pmm = 
(1/2 — /3 + S) of the direction dependent term. At present the most stringent 
limit on that parameter is \Pmm \ < 3.2 x 10~ 9 [6, 18] determined over 23 years 
ago in an outstanding experiment whose uncertainty has not been approached 
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since. Our experiment is the first one that can confirm that result with roughly 
equivalent uncertainty (4.2 x 10~ 9 ). Kcnncdy-Thorndike experiments [7, 8, 9] 
measure the coefficient Prt = ((3 — a — 1) of the velocity dependent term. 
The most stringent limit on \Pkt\ has been recently improved from [8] by a 
factor 3 to \Pkt\ < 2.1 x 10~ 5 [9]. We improve this result by a factor of 30 to 
\Prt\ < 6.9 x 10 -7 . Finally clock comparison and Dopplcr experiments measure 
a, currently limiting it to \a + 1/2] < 8 x 1CT 7 [19, 20, 21]. The three types 
of experiments taken together then completely characterize any deviation from 
Lorentz invariance in this particular test theory. 

Our cryogenic oscillator consists of a sapphire crystal of cylindrical shape 
operating in a whispering gallery mode (see fig. 1 for a schematic drawing and 
[22, 23] for a detailed description). Its coordinate frequency can be expressed 
by v = m/t c where t c is the coordinate travel time of a light signal around the 
circumference of the cylinder (of radius r) and m is a constant. From (2) the 
relative frequency difference between the sapphire oscillator and the hydrogen 
maser (which realizes coordinate time in S) is 

^ = p kt v -^ + Pmm V ^ r^e { t^ + o( c -3) (3) 

Vq c c Ztt J q 

where v$ = m/(27rr/c), v (t) is the (time dependent) speed of the lab in S, and 
tp is the azimuthal angle of the light signal in the plane of the cylinder. The 
periodic time dependence of v and 9 due to the rotation and orbital motion 
of the Earth with respect to the CMB frame allow us to set limits on the two 
parameters in (3) by adjusting the periodic terms of appropriate frequency and 
phase (see [24] for calculations of similar effects for several types of oscillator 
modes). Given the limited durations of our data sets (< 9.5 days) the dominant 
periodic terms arise from the Earth's rotation, so retaining only those we have 
v(i) = u + fi x R with u the velocity of the solar system with respect to the 
CMB, fi the angular velocity of the Earth, and R the geocentric position of the 
lab. We then find after some calculation 

A^/^o = PKT(HsinX) 

+Pmm (AiosA + i?cos(2A) + CsinA + DsinAcosA + _EsinAcos(2A)) 

(4) 

where A = Qt + <j), and A-E and (j) are constants depending on the latitude 
and longitude of the lab (w 48.7°N and 2.33°E for Paris). Numerically H w 
-2.6 x 10" 9 , A w -8.8 x 10~ 8 , B w 1.8 x 10~ 7 , C-E of order 10" 9 . We note 
that in (4) the dominant time variations of the two combinations of parameters 
are in quadrature and at twice the frequency which indicates that they should 
decorelate well in the data analysis allowing a simultaneous determination of 
the two (as confirmed by the correlation coefficients below). Adjusting this 
simplified model to our data we obtain results that differ by less than 10% from 
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the results presented below that were obtained using the complete model ((3) 
including the orbital motion of the Earth) . 

The cryogenic sapphire oscillator (CSO) is compared to a commercial (Da- 
tum Inc.) active hydrogen maser whose frequency is also regularly compared 
to caesium and rubidium atomic fountain clocks in the laboratory [25]. Both 
oscillators are operated in temperature controlled rooms, with the temperature 
sensitive electronics mounted on an actively temperature stabilized panel. The 
CSO resonant frequency at 11.932 GHz is compared to the 100 MHz output 
of the hydrogen maser. The maser signal is multiplied up to 12 GHz of which 
the CSO signal is subtracted. The remaining w 67 MHz signal is mixed to a 
synthesizer signal at the same frequency and the low frequency beat at w 64 Hz 
is counted, giving access to the frequency difference between the maser and the 
CSO. The instability of the comparison chain has been measured and does not 
exceed a few parts in 10 16 . The typical stability of the measured CSO - maser 
frequency after removal of a linear frequency drift is shown in fig. 1. 

Our experimental data consists of seven sets of measured values of Av/vq 
of varying length (2.8 to 9.5 days, 37 days in total) taken in Nov. /Dec. 2001 
and Mar. /Apr. /Sep. 2002 (see fig. 2). The sampling times to are generally 
100 s except for the first two sets for which to = 12 s and 5000 s respectively. 
To analyze our data we simultaneously adjust an offset and a rate (natural 
frequency drift, typically «2x 10~ 18 s _1 ) per data set and the two parameters 
of the model (3), a total of 16 parameters. In the model (3) we take into account 
the rotation of the Earth and the Earth's orbital motion, the latter contributing 
little as any constant or linear terms over the durations of the individual data 
sets are absorbed by the adjusted offsets and rates. To ensure homogeneity 
between the data sets we average all sets to t = 5000 s. 

We first carry out an ordinary least squares (OLS) adjustment obtaining 
\Pmm\ = (-2.8 ±3.6) x 10~ 9 and \P KT \ = (-4.3 ±2.8) x 10~ 7 . The correlation 
coefficient between the two parameters is less than 0.01 indicating that the two 
are indeed well decorelated and can be determined simultaneously. All other 
correlation coefficients between either of the two parameters and the 14 adjusted 
offsets and rates are less than 0.07. We note, however, that the residuals have 
a significantly non-white behavior as one would expect from the slope of the 
Allan deviation of fig. 1. The power spectral density (PSD) of the residuals 
when fitted with a power law model of the form S y (f) — kf^ yields typically 
fi « —1.5. In the presence of non- white noise OLS is not the optimal regression 
method [26, 27] as it can lead to significant underestimation of the parameter 
uncertainties [26]. 

An alternative method is weighted least squares (WLS) [27] which allows 
one to account for non-random noise processes in the original data by pre- 
multiplying both sides of the design equation (our equation (3) plus the 14 
offsets and rates) by a weighting matrix containing off diagonal elements. To 
determine these off diagonal terms we first carry out OLS and adjust the S y (f) = 
kf^ power law model to the PSD of the post-fit residuals determining a value 
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of \i for each data set. We then use these \i values to construct a weighting 
matrix following the method of fractional differencing described, for example, 
in [26]. The WLS regression yields \Pmm\ = (1.5 ± 3.1) x 1(T 9 and \P KT \ = 
(—3.1 ± 3.7) x 10~ 7 (lcr uncertainties), with the correlation coefficient between 
the two parameters less than 0.01 and all other correlation coefficients < 0.07. 
The best fit power law of the WLS residuals is now compatible with /i ss 0. The 
inset of fig. 2 shows the last data set (Sep. 2002) and the best fit model after 
pre-multiplication of both by the weighting matrix. 

The two methods give similar results for the two parameters but we consider 
WLS as more reliable and take its outcome as the final results of our statistical 
analysis. 

Systematic effects at diurnal or semi-diurnal frequencies with the appropri- 
ate phase could partially cancel a putative sidereal signal as our total data span 
(« 293 days) allows only partial separation of the diurnal from the sidereal sig- 
nal. The statistical uncertainties of Pmm and Pkt obtained from the WLS fit 
above correspond to sidereal and semi-sidereal terms (from (4)) of « 1 x f0~ 15 
and ~6x f0~ 16 respectively so any systematic effects exceeding these limits 
need to be taken into account in the final uncertainty. We expect the main 
contributions to such effects to arise from temperature, pressure and magnetic 
field variations that would affect the hydrogen maser, the CSO and the associ- 
ated electronics, and from tilt variations of the CSO which are known to affect 
its frequency. Measurements of the tilt variations of the CSO show amplitudes 
of 4.6 /Ltrad and 1.6 /zrad at diurnal and semi-diurnal frequencies. To estimate 
the tilt sensitivity we have intentionally tilted the oscillator by w 5 mrad off its 
average position which led to relative frequency variations of ks 3 x 10~ 13 from 
which we deduce a tilt sensitivity of « 6 x 10 _17 /xrad _1 . This value corresponds 
to a worst case scenario as we expect a quadratic rather than linear frequency 
variation for small tilts around the vertical. Even with this pessimistic estimate 
diurnal and semi-diurnal frequency variations due to tilt do not exceed 3 x 10~ 16 
and f x fO -16 respectively and are therefore negligible with respect to the sta- 
tistical uncertainties. The temperature sensitive electronics were mounted on 
an actively temperature controlled panel reducing temperature fluctuations by 
about one order of magnitude. Temperature measurements of the CSO lab and 
the electronics panel taken during some of the experimental runs show room 
temperature variations with amplitudes of 0.3 °C and 0.1 °C for the diurnal 
and semi-diurnal components which are reduced to 0.04 °C and 0.01 °C on 
the panel. The hydrogen maser is kept in a dedicated clock room with tem- 
perature variations below the above values. Measurements of magnetic field, 
temperature and atmospheric pressure in that room and the maser sensitivi- 
ties as specified by the manufacturer allow us to exclude any systematic effects 
on the maser frequency that would exceed the statistical uncertainties above. 
Switching off the temperature stabilization of the electronics panel shows no 
discernible effect so we are confident in excluding any systematic effects from 
that source. When heating and cooling the CSO lab by w 3°C we see frequency 
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variations of w 5 x 10~ 15 per °C. From the temperature measurements dur- 
ing the experimental runs we therefore deduce a total diurnal and semi-diurnal 
effect of « 1.5 x 1CP 15 and « 5.0 x 1CP 16 respectively. We assume the pes- 
simistic scenario where all of the diurnal and semi-diurnal effect is present at 
the neighboring sidereal and semi-sidereal frequencies which leads (from (3)) to 
uncertainties from systematic effects of ±5.8 x 10~ 7 on Pkt and ±2.8 x 10 -9 
on Pmm ■ We note that the phase of the perturbing systematic signal will vary 
over the course of our measurements due to natural causes (meteorology, day- 
time changes etc.) and to the sidereal/diurnal frequency difference so our final 
uncertainties given below are the quadratic sums of the above values and the 
statistical uncertainties from the WLS adjustment. 

In summary, we have reported an experimental test of Lorentz invariance 
that simultaneously constrains two combinations of the three parameters of the 
Mansouri and Sexl test theory (previously measured individually by Michelson- 
Morley and Kenncdy-Thorndikc experiments). Our experiment limits \S — (3 + 
1/2 1 < 4.2 x 10~ 9 which for the first time confirms the 1979 value of [6], and 
\(3 — a — 1| < 6.9 x 10~ 7 which improves the best present limit [9] by a factor of 
30. As a result, the Lorentz transformations are confirmed in this particular test 
theory with an overall uncertainty of < 8 x 10~ 7 limited now by the determina- 
tion of a from Dopplcr and clock comparison experiments [19, 20, 21]. This is 
likely to be improved in the coming years by experiments such as ACES (Atomic 
Clock Ensemble in Space [28]) that will compare ground clocks to clocks on the 
international space station aiming at a 10 fold improvement on the determina- 
tion of a. We hope to improve our experiment by a better characterization and 
control of systematic effects (in particular due to temperature) and by further 
data integration. This should allow us to improve our limits by another factor 
2 or 3 in the near future. 
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Fig. 1. Typical relative frequency stability of the CSO - H-maser difference 
after removal of a linear frequency drift. The inset is a schematic drawing of 
the cylindrical sapphire oscillator with the Poynting vector P in the whispering 
gallery (WG) mode, the velocity v(£) of the cylinder with respect to the CMB, 
and the relevant angles for a photon in the WG mode. 
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Fig. 2. Measured relative frequency difference between the CSO and the hydro- 
gen maser after removal of the adjusted offsets and rates, and best fit model (3) 
of Lorentz invariance violation (solid line). The inset shows the last data set 
after pre-multiplication of the data and the model by the weighting matrix of 
the WLS adjustment (for presentation purposes data points were averaged by 
a factor 2). The error bar is the maximum signal amplitude within the final la 
uncertainties of the parameters. 
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